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We introduce a model describing the paths that pin an elastic interface moving in a disordered 
medium. We find that the scaling properties of these "elastic pinning paths" (EPP) are different 
from paths embedded on a directed percolation cluster, which are known to pin the interface of the 
"directed percolation depinning" class of surface growth models. The EPP are characterized by a 
roughness exponent a = 1.25, intermediate between that of the free inertial process (a = 3/2) and 
the diode-resistor problem on a Cayley tree (a — 1). We also calculate numerically the mean cluster 
size and the cluster size distribution for the EPP. 
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The problem of interface roughening in the presence of 
quenched disorder is a topic of recent interest, due to its 
importance as a paradigm in condensed matter physics 
and due to the broad range of applications [fjj . In a typ- 
ical case, the interface moves in a (d + l)-dimensional 
disordered medium driven by a homogeneous force F. 
At small forces, the interface is pinned by the impurities 
of the medium, while the interface undergoes a depin- 
ning transition at a critical force F c , and for F > F c the 
interface moves with a nonzero velocity. The spatial fluc- 
tuations of the interface are characterized by the scaling 
of the saturated interface width W sa t (L) with the system 
size, 



W sat (L) ~L C 



(1) 



where a is the roughness exponent. 

It has been proposed EqJ3] that the depinning transi- 
tion can be described by the following equation of motion 
for the interface height y(x, t) 



d_ 

at 



y{x, t) = V 2 y + r)(x, y) + F, 



(2) 



where x is the d-dimensional coordinate parallel to the 
interface. The first term on the right hand side of (||) 
represents the surface tension favoring a smooth inter- 
face, and we say that the interface is elastic. The second 
term is a random field that mimics the quenched disor- 
der of the medium and is assumed to have zero mean and 
short-range correlations. 

The universality class corresponding to Eq. (Q)is called 
quenched Edwards- Wilkinson (QEW), because (||) is sim- 
ilar to the Edwards- Wilkinson equation The differ- 
ence, which changes the behavior of (g) drastically, is the 
presence of spatially dependent quenched disorder rj{x, y) 
instead of time-dependent shot noise rj(x,t). Numerical 



studies J5|-[7j of the depinning transition yield a rough- 
ness exponent a ~ 1.25 in d = 1 ((1 + 1) dimensions) 
and a ~ 0.75 for d = 2 ((2 + 1) dimensions). These val- 
ues are lower than the results of perturbation theory M , 
which yields a = 3/2 in d = 1 and a = 1 for d = 2. 
On the other hand, the numerical values are significantly 
higher than the prediction of a functional renormaliza- 
tion group treatment which gives a ~ 1 and a ~ 2/3 for 
d = 1 and d = 2, respectively ||. 

The relevance of directed percolation to interface de- 
pinning has been established for a different class of 
models, called directed percolation depinning models |)), 
which are in the same universality class as Eq. (g) when a 
Kardar-Parisi- Zhang (KPZ) @ term A(Vy) 2 is included. 
In these models, the interface is pinned by paths on a di- 
rected percolation cluster of pinning sites ||. Thus, the 
scaling properties of the interface at the depinning tran- 
sition in (1 + 1) dimensions can be obtained by a map- 
ping onto directed percolation (DP) [|To|JlT[l , from which 
a roughness exponent a ~ 0.63 is obtained. 

In this paper, we consider the paths which pin the in- 
terface for the QEW universality class for the case d = 1. 
We term these paths elastic pinning paths (EPP). We ap- 
ply the concepts of directed percolation to investigate the 
scaling properties of the EPP. Our numerical results pro- 
vide an independent check for the anomalous roughness 
exponent a ~ 1.25 obtained with the models of the QEW 
universality class. We also consider two known random 
walk models, which yield a = 3/2 and a = 1, respec- 
tively. By comparing the EPP to these random walks we 
obtain some insight why the roughness exponent of Eq. 
(g) lies in the interval 1 < a < 3/2. 

To motivate the definition of the EPP, we first con- 
sider a discrete solid-on-solid model corresponding to Eq. 
(g) The interface position hi is defined on a square 
lattice of lateral size L. We assign to each site on the 
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lattice a random number rji t h which can have two val- 
ues, r)i t h = 1 (unblocked cell) with probability p, and 
r]i t h = —1 (blocked cell) with probability 1 — p. A local 
force is defined by 

fi = hi+i + hi-i - 2hi + r) i:hi . (3) 

At time t = the interface is flat, and at a given time 
the height of the i-th column is increased by one if the 
local force fi is positive. 

At a critical value of the probability p = p c , the inter- 
face is pinned by one of the pinning paths. According to 
the dynamical rules, a pinned interface satisfies fi < 
for all i. We define the increments of a given path as 
Ai = hi — hi-i, so according to Eq. (||) a spanning path 
stopping the QEW interface satisfies 

A l+ i < A, - m th . (4) 

By induction, one can show that in addition to the con- 
dition (H), a lower bound for A J+ i — Ai holds at every 
time step of the interface evolution 

A 4+ i - A, > -2. (5) 

Equations (0) and (Q) define the possible pinning 
paths. Note that paths in the same cell (i, h) can have 
different increments Aj. We start at i = 1 with hi = 
and initial increment Ai — 0. The paths in column i + 1 
are updated according to the following three rules: 

Rule (i): If the cell (i,h) is blocked (r)i t h = — 1), 
then the path is splitted into four paths, where the po- 
sitions at i + 1 are hi+i = hi + Aj + 1, hi+\ = hi + Ai, 
hi+i =hi+Ai-l, and h i+1 = h { + A. t - 2 

Rule (ii): If the cell (i,h) is unblocked, rj^h = l,we 
have: /ij+i = hi + Aj — 1, and /ij+i — hi + Ai — 2. 

Rule (iii): The path stops when hi < 0. 

After moving to the new cell the increment Aj+i is 
updated and the rules are applied again. Rules (i) and 
(ii) are the implemetation of Eq. (||). Rule (iii) is moti- 
vated by the fact that if the path deviates too much in 
the downward direction, it would not have a chance to 
block the growth since, in a system with periodic bound- 
ary conditions, the path should return to the same point 
where it starts. In Fig. [I] we show a typical set ("clus- 
ter" ) of directed paths. The paths are characterized by 
large local slopes which is the main feature of the QEW 
interface at the depinning transition. 

The scaling properties of the directed paths require two 
characteristic lengths, £|| and the correlation length 
parallel to and perpendicular to the preferred direction 
of the paths pl| . The correlation lengths diverge at the 
critical concentration of pinning centers p c as 

£n(p)~b-p c n u(p)^\p-po\- ux , (6) 

where vu and is± are two different universal exponents 
due to the anisotropy given by the preferred direction of 
the paths. 



A cluster consisting of pinning paths is defined by all 
paths generated by the Rules (i— iii) for one realization of 
the disorder. Let us denote by s the number of sites in a 
given cluster. The cluster size distribution n s {p) 1 defined 
as the average number of clusters of s sites per lattice 
site, shows a power-law behavior at p c . For p < p c only 
finite clusters are present, so there exists an effective cut- 
off for the cluster size, sq ~ \p — p c \~ x l a ■ The cluster size 
distribution has the scaling form n s (p) ~ s~ T g(s/so), 
where g{x) is a scaling function that decreases faster 
than a power-law for x 3> 1. The mean cluster size (s) 
also diverges at p c as a power-law (s) ~ \p — p c \ f , with 
7 =(2-T)/<r. 

The scaling relations presented so far are valid for in- 
finite lattices. Finite size scaling considerations allow us 
to write § W(L) ~ £± ~ ^ x/l/ " ~ L"^ V K Hence from 

a = v±/v\\. (7) 

In our simulations we compute all the exponents char- 
acterizing the scaling behavior of the EPP. Results are 
shown in Fig. [|. We calculate the correlation length ex- 
ponents and find v\\ ~ 1.33 and v±_ ~ 1.67 (see Table |). 
Using (Q) we find a ~ 1.25 in agreement with the expo- 
nent found for the QEW interface . The study of the 
mean cluster size yields an exponent 7 ~ 2.43, while the 
exponent of the cluster size distribution is r ~ 1.43. 

Next we calculate the scaling of the mean square fluc- 
tuations or "width" of the paths defined as {hf) 1 ^ 2 as a 
function of the parallel coordinate i calculated at p c . We 
find (h 2 ) 1 ' 2 ~ i a with a ~ 1.27, in agreement with the 
exponent a = v_\_/v\\ = 1.25 that we find using (j7j) and 
our numerical results for v±_ and v\\ . 

We also compute the mean square fluctuations of the 
increments Aj. We find (A 2 ) 1 / 2 ~ i a where a' = a — 1 
because hi is the integral of Aj along the path, so by 
integration the exponent a is 1 + a' . 

It is interesting to note that Rules (i) and (ii) can 
be modified without changing the universality class. In- 
deed, instead of four choices for Aj+i in Rule (i) and 
two choices for Aj+i in (ii), we can set Aj+i = Aj + k, 
k = -1,0,1 (Rule (ia)), and A l+ i = A 4 -l, (Rule (iia)), 
or even Aj+i = A^ + fc, with k = — 1, 1 (Rule (ib)). These 
changes are analogous to the change of coordination num- 
ber in directed percolation which change the value of the 
critical threshold but do not change the universality class. 
The results in Table || correspond to the most extensive 
simulations which are performed with the Rules (ib) and 
(iia). 

Next, we discuss the relation of the EPP to two known 
universality classes — which can be considered as random 
walk models — and to which the EPP can be modified by 
changing the interaction with the disorder. For the first 
universality class we assume that the noise is determined 
by the position (i, Ai) in the increment space, instead 
by the real space position (i, hi) for the EPP. Since the 
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average span of Aj ~ i a is much smaller than the aver- 
age span of hi ~ i a , many of the EPP that were treated 
differently (since they have different hi coordinate) be- 
come indistinguishable, since many of them can have the 
same value of Aj, so that the number of different paths 
decreases significantly compared to the EPP. 

The span of the cluster is determined by the top-most 
trajectory in real space which is in turn represented by 
the top-most trajectory in the increment space (i,Aj). 
Since the noise is chosen according to the position in 
the (i,Aj) space, two paths which arrive to the same 
point (zojAo) remain together because also the noise is 
the same. Thus, two given paths cannot cross each other 
in the plane (i, Aj) (see Fig. ||) and the top-most trajec- 
tory in the increment space is composed by a unique and 
well-defined path. This path is a simple random walk in 
the {i, A,) plane, which has a certain probability of go- 
ing up and down depending on the noise. Therefore the 
average span scales as Aj ~ i a with a' — 1/2. With 
q = a' + 1 we get a = 3/2 > «epp = 1.25. This is in 
agreement with the fact that the standard deviation of 
the free inertial process |l4|,[l5| increases as i 3 ^ 2 . 

The critical probability, p c , now corresponds to the un- 
biased random walk, and can be readily computed ana- 
lytically for each of the modification of the Rules (i) and 
(ii). The critical probabilities are always larger that the 
corresponding values in the EPP case, being p c = 1/2 for 
the rules (ib) and (iia). The exponent t now is related 
to the probability of the first return to the origin after 
some number of steps i. For the free inertial process this 
probability decays as i -5 / 4 which differs from the 
random walk result i -3 / 2 . This indicates that the top- 
most trajectory in the (z, Aj) plane determines only the 
scaling of the span of the clusters, but other properties, 
such as the probability of first return to the origin, are 
not determined by the top-most trajectory. 

For the EPP, where the noise is chosen from the po- 
sition in real space (i,h), the top- most trajectory in the 
(i, Aj) plane still determines the scaling properties of the 
span on the clusters. However, in this case, the top-most 
trajectory in the increment space can be composed by 
parts of several different paths (Fig. ||). This is so be- 
cause two paths can cross each other in the (i, Aj) plane 
because the noise is determined by the position in the 
(i, h) plane which can be different for the two paths at 
the crossing point (io, A ) (Fig. ||). This is the main dif- 
ference between the EPP and the free inertial process. 

The second universality class to which the EPP reduces 
is obtained by defining the value of the noise for each path 
independently. Even for two paths that meet at a point 
(i,hi) with the same Aj, the values of the noise are cho- 
sen as independent of each other. Since now all paths are 
independent and there are no loops, this model can be 
exactly mapped — with respect to Aj and i — to the diode 
resistor Cayley tree problem solved in |L7| . The i coordi- 
nate is identified with the time coordinate of the Cayley 
tree cluster which, in turn, corresponds to the chemical 
or minimum path in the longitudinal hyperplane of the 



Cayley tree cluster (see Q), and has v\\ = 1/2 @. The 
coordinate Aj is the remaining transversal coordinate of 
the Cayley tree cluster and has uj_ = p7| . Then, we find 
a' = 0, and a = a' + l = l< olepp = 1-25. We also find 
numerically a' — 0. The critical probability p c in this 
case can be computed analytically, and it decreases in 
comparison with the EPP case. For rules (ib) and (iia) 
Pc = 1/4. 

What can we learn from the comparison of the EPP to 
the two discussed random walk models? Since the EPP 
can cross each other in the (i, Aj) plane, there is no 
unique top-most path but a top-most trajectory consists 
of several paths. By being at some point not the top- 
most path, a given path can "optimize" its way through 
the randomness in the sense that it visits more blocked 
cells. Thus, it can stay alive longer (hi > 0) compared 
to the free inertial process where the unique top-most 
path determines the scaling properties. This means that 
the EPP has a smaller roughness exponent than the free 
inertial value of 3/2. 

The other random walk model, where the noise of dif- 
ferent paths is always independent, has a roughness ex- 
ponent a = 1 which is even lower than that of the EPP. 
In order to understand the reason, we note that the paths 
have at each point many more possibilities than the EPP. 
Thus, it is not surprising that there are always paths 
which stay alive even for large i. It seems reasonable 
that the scaling properties are dominated by these paths, 
causing a smaller roughness exponent. 

To summarize, we introduce the EPP model to de- 
scribe the paths that pin an elastic interface moving in 
a disordered medium, and we find that it is in a uni- 
versality class different from directed percolation paths 
that are responsible for pinning of the class of models 
of the directed percolation depinning universality class. 
The critical paths that are characterized by an exponent 
olepp — 1.25, describe the scaling properties of an elastic 
interface pinned by the quenched disorder of the medium. 
The roughness exponent of the EPP lies between that of 
the free inertial process with a = 1.5 juj, and that of 
the diode resistor Cayley tree problem with a = 1 [ p7[ . 
The comparison of the EPP with these two random walk 
models sheds some light on the value of a = 1.25. 
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TABLE I. Critical exponents for the EPP presented in 
this paper along with the results for the elastic interface in 
the QEW universality class, and the DP universality class. 
We find Pc = 0.282 ± 0.001 for the EPP. 





EPP 


a 


1.26 ±0.03 


v \\ 


1.33 ±0.04 


v± 


1.67 ±0.04 


T 


1.43 ±0.02 


7 


2.43 ±0.05 



QEW (interface) 



1.25 ±0.01 
1.35 ±0.04 
1.68 ± 0.04 



DP 



0.63 ± 0.04 M 
1.73 ±0.01 m\ 
1.09 ±0.01 L] 
1.28 ±0.02 [L] 
2.28 ±0.01 y] 




FIG. 1. Example of one cluster of pinning paths gener- 
ated by the proposed rules. The sites in a square lattice are 
unblocked (n — 1, shown in white) with probability p and 
blocked (n — —1, shown as cross-hatched) with probability 
1 — p. The path starts at xi = with hi = and Ai = and 
we apply the rules (ib), (iia), and (iii). 
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FIG. 2. Log-log plots of the mean cluster size < s > 
(shifted for clarity), the correlation lengths £x and £n in the 
perpendicular and parallel directions as a function of the re- 
duced probability (p c —p)/p c . Simulations are averaged over 
10 7 clusters. 



(a) EPP model 



(b) Free inertial model 
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FIG. 3. Schematic illustration of the top-most trajectory 
in the (i, Ai) plane for (a) the EPP, and (6) the modification 
of the EPP where the noise is chosen according to the position 
of the path in the (i, Ai) plane (free inertial process). In both 
figures we plot two paths (solid and dashed lines) which inter- 
sect at the point (io, Ao). In the EPP case the paths can cross 
each other, because the noise is determined by the position 
in real space (i,h). Thus, also paths in real space can cross 
each other. In the free inertial process, however, if two paths 
intersect at (io,Ao) then they continue together. Therefore, 
paths in real space cannot cross each other. 
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